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1.0 Introduction : 

CLFE2D was developed under the assumption of a state of generalized 
plane strain. The most general displacement field for such problems can 
be written 


“total * u <y> 2 ) + * 

(la) 

V * V(y,z) 

(lb) 

w = W(y,z) 

(lc) 


The assumption of generalized plane strain allows the user to study the 
response of an infinitely long body of arbitrary cross-section to 
various loadings with the restriction that all quantities except axial 
displacement are independent of the axial coordinate x* The coordinate 
system used by CLFE2D is shown in Fig, (1,1), The x, y and z axes 
represent the global coordinate system while the material principal 
directions are represented by the 1, 2 and 3 axes. The region of 
interest to be modeled by finite elements is then a typical cross 
section (or a symmetric portion of it). 

The element used in CLFE2D is a four node, linear, general quadri- 
lateral, i soparametric element. A three node version of the element can 
be obtained by “double-numbering" the first local node in the element 
connectivity data. The available element geometries are shown in 
Fig. (1.2). CLFE2D requires that the local node numbering scheme, for a 
given element, proceed in a counter-clockwise direction. 

An in-depth description of the finite element used in this program 
can be found in reference [1]. The finite element formulation for the 
generalized plane strain problem can be found in reference [2], 

In addition to the above overview of the program capabilities, this 
report presents the input requi rements , the output capabilities, details 





of the Input card sequence, a note on symmetry requirements and an 
example problem. 

2.0 CLFE2P Input Requirements : 

In this section a description of each type of card image in the 
input file is given. The required format of each ctrd image is given in 
section 3.0* 

Type 1 Title (Cards 1 & 2) 

The title consists of two 80 character lines which will be 
printed on the output file. 

Type 2 Master control Information (Card 3) 

On this card the user specifies the number of elements and 
nodes in the mesh, the number of different angles of material 
principal orientation, the number of specified nodal degrees of 
freedom and forces. The user also specifies the desired plot 
option and a data check option on this card (see section 

3.0) . The data check option permits the user to check the 
input vithout proceeding to the solution. 

Type 3 Nodal Coordinate Scale Factors (Card 4) 

To ease input the user can specify scale factors which 
independently scale the magnitude of the Y and Z nodal 
coordinate values. 

Type 4 Printer Control Keys (Card 5) 

On this card the user sets the various output keys (see section 

3.0) which allow for the generation of a tailored output file. 
Type 5 Plotter Control Information (Card 6) 

On this card the user specifies the plotter scale factors (see 


6 




Output Capabilities). Card 6 can be omitted if output plots are 
not desired. 

Type 6 Material Property Information (Card 7 and 8} 

On these two cards the user inputs one set of three dimensional 
orthotropic material properties. These properties must include 
thermal or hygroscopic expansion coefficients if the response 
to hygrothermal loading is desired. 

Type 7 Material Property Orientation Information (Card 9) 

The user inputs an angle number and magnitude (degrees) for 
each different material property orientation in the laminate. 

Type 8 Element Data Information (Card 10) 

The element connectivity data, material property orientation 
angle number, and the output points (see section 3.0) for each 
element are specified on this card. Card 10 is repeated NEM 
times, where NEM is the number of elements in the mesh. 

Type 9 Nodal Data Information (Card 11) 

The nodal coordinates are specified on this card. The nodal 
coordinate values must be greater than or equal to zero. Card 
11 is repeated NODS times, where NODS is the number of nodes in 
the mesh. 

Type 10 Specified Nodal Displacement Information (Card 12) 

On this card the user identifies the node number and displace- 
ment of each constrained degree of freedom. Card 12 is 
repeated for each degree of freedom that has a specified dis- 
placement. If there are no constrained degrees of freedom, 

Card 12 is omitted from the input file. 


27 
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Type 11 Specified Nodal Force Information (Card 13) 

On this card the user identifies the node number and the compo- 
nent magnitude of each of the specified nodal forces. Card 13 
is repeated for each component of the specified nodal forces. 

If there are no specified nodal forces, Card 13 is omitted from 
the input file. 

Type 12 Normal Strain and Hygrothermal Load Information (Card 14) 

Here the user specifies the applied normal strain e° along with 

a 

the temperature or moisture change. It is possible to specify 
these values as zero. 

3.0 CL FE2P Output Capabilities : 

CLFE2D solves for the nodal displacements, strains, stresses and 
strain energy densities for each element in the mesh. These values are 
calculated in the global (x, y, z) coordinate system. For a given 
element, the values of strain, stress and strain energy density are 
calculated at a combination of any nine output points. These output 
points are specified for each element on Card 10 of the input file. 
Figure (1.3) illustrates the location of these output points on a given 
element. Output point 5 is the only valid output point for a 
triangular element. Specifying other values for triangular elements can 
cause errors. 

It should be noted that the strains, stresses, and strain energies 
that are printed are the response to the combined mechanical and 
Hygrothermal loadings. The u displacement that is printed is the U(y,z) 
of eqn. (la). To obtain the total u displacement at a desired x 
location, one must take into account the applied normal strain e°. 

A 
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In CIFE2D the user has the ability to tailor the output file. By 
setting various printer control keys on Card 5 of the input file, he can 
opt to have all or any of the following printed; 


Key (1) 

Element Data 

Key (2) 

Nodal Data 

Key (3) 

Specified Nodal Displacements 

Key (4) 

Specified Nodal Forces 

Key (5) 

Displacements 

Key (6) 

Strains, Stresses and Strain Energy Densities 


The user can also opt to plot the undeformed and deformed meshes by 
selecting the appropriate plot option on Card 3. The user controls the 
scale of these plots by specifying various plot scale factors on Card 
6. The various plot options and scale factors are described in detail 
in the following outline. 

1. Plot Option NPLOT (Card 3) 

NPLOT is one of the variables which is user specified on Card 3. 

The user determines which plots will be printed by setting NPLOT to 
any of the following values. 

0 = No plots 

1 = Undeformed plot only 

2 = Deformed plot only 

3 * Both undeformed and deformed plots 

2. Plotter Control Information (Card 6) 

On Card 6 of the input file, the user specifies the following plot 
scale factors: 


PYSCL * Plot Y-direction scale factor 

This varlaole determines the maximum length of the plot In the 
Y-direction. PYSCL must be chosen such that the maximum Y- 
coordinate In the mesh will be scaled to a value less than or 
equal to 18 inches. 

PZSCL = Plot Z -direction scale factor 

This variable determines the maximum length of the plot in the 
Z -direction. PZSCL must be chosen such that the maximum Z- 
coordinate in the mesh will be scaled to a value less than or 
equal to 9 inches. 

VMAX « Maximum allowable v-d^splacement 

This variable determines the size of the maximum v-displacement 
in inches. All of the smaller v-displacements are scaled 
proportionally to this value, 

WMAX « Maximum allowable w-displacement 

This variable determines the size of the maximum w-displacement 
in inches. All of the smaller w-displ acements are scaled 
proportionally to this value. 
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4.0 CLFE2D Input Card Sequence 


Cards 1 & 2 (20A4) 

Column 

1-80 


Title Cards 
Contents 


Title 


Card 3 (1615) 
““Column 


1-5 

NEM 

6-10 

NODS 

11-15 

NANG 

16-20 

NSDF 

21-25 

NSBF 

26-30 

NPLOT 


41-45 


Master Control Card 
Contents 

Number of El ements 
Number of Node Points 
Number of Different Angles 
Number of Specified Nodal Displacements 
Number of Specified Nodal Forces 
Plot Option 

* 0 none 

= 1 undeformed plot only 
= 2 deformed plot only 

* 3 both undeformed and deformed 
NCHECK Data Check Option 


By setting NPLOT .EQ.l and NCHECK .NE.O a data check and an 
undeformed plot will be generated. To obtain the complete 
solution, another run must be made with NCHECK .EQ.Q. 


Card 4 (8D10.5) 

Scale Factors 

Column 


Contents 

1-10 

SCAY 

Y- Nodal Coordinate Scale Factor 

11-20 

SCAZ 

Z-Nodal Coordinate Scale Factor 

Card 5 (1615) 


Printer Control Card 

Col umn 


Contents 

1-5 

KEY ( 1 ) 

Key for Printing Element Data 

6-10 

KEY(2) 

Key for Printing Nodal Data 

11-15 

KEY(3) 

Key for Printing Specified Nodal 



Displacements 

16-20 

KEY(4) 

Key for Printing Specified Nodal Forces 

21-25 

KEY(5) 

Key for Printing Displacements 

26-30 

KE Y ( 6 ) 

Key for Printing Strains, Stresses and 


Strain 

Energies Per Unit Volume 

(if 

KEY(i) .NE 

'. 0 - Print) 

(NPLOT .EQ. 0 

skip Card 6 and go to Card 7) 

Card 6 (8E10.5) 

Plotter Control Card 

Column 


Contents 

1-10 

PYSCL 

PI ot Y-Scale Factor 

11-20 

PZSCL 

Plot Z-Scale Factor 

21-30 

VMAX 

Maximum Allowable V -Displacement 

31-40 

WMAX 

Maximum Allowable W-Displacement 
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Card 7 (8D10.5) 

“Column 

1-10 HOP(l) 

11-20 PROP ( 2 ) 

21-30 PROP ( 3 ) 

31-40 PR0P(4) 

41-50 PROP ( 5 ) 

51-50 PROP ( 6 ) 

Card 8 (8D10.5) 

Column 

f^lO~ Prop(7) 

11-20 Prop(8 

21-30 Prop(9) 

31-40 Prop 10) 

41-50 Prop (ll) 

51-60 Prop(12) 


Material Properties Card A 
Contents 
E-TI 
E-22 
E-33 
G-23 
6-13 
6-12 

Material Properties Card B 
Contents 
NU-23 
NU-13 
NU-12 
Alpha-11 
Alpha-22 
A1 pha-33 


Card 9 (8D10.5) 
Column 

l'-lO Ang(l) 

11-20 Ang(2) 


Angle Data Card 
Contents 

Angle No. 1 (in Degrees) 
Angle No. 2 (in Degrees) 




Ang(NANG) 


Angle No. NANG (in Degrees) 


Card 10 (5X.6I5) 
Column 


Element Data Card(s) 
Contents 


1-5 


6-10 

11-15 

16-20 

21-25 

26-30 

31-35 


Blank (Element Numbers may be inserted for 

Reference) 

NOD ( N , 1 ) Node #1 of Element N 

N0D(N,2) Node # 2 of Element N 

N0D(N,3) Node #3 of Element N 

N0D(N,4) Node H of Element N 

IANG(N) Angle Number of Element N 

I$TRS(N,1) Stress and Strain Output Code for Element N. 


■ 0 - None 
= 1 - Side 1 
= 2 - Side 2 
= 3 - Side 3 
= 4 - Side 4 
- 5 - Center 



= 6 

- "Node 

1" 


= 7 

- "Node 

2" 


= 8 

- "Node 

3" 


= 9 

- "Node 

4" 

36-40 

ISTRS ( 

;n, 2 ) 


41-45 

ISTRS( 

!n,3) 



71-75 ISTRS ( N , 9 ) 

*Repeat Card 10 NEM Times 


Card 

11 (5X,2D10.5) 

Nodal Data Card(s) 


Column 


Contents 


1-5 

Blank 

(Node Numbers may be inserted for Reference) 


6-15 

Y(N) 

Y-Coordinate of Node N (This value must 




be > zero) 


16-25 

Z(N) 

Z -Coordinate of Node N (This value must 




be > zero) 



^Repeat Card 11 NODS Times # 



# If NSDF 

.EQ. 0 - Skip Card 12 and Go on to Card 13 

Card 

12 (215,010. 5) 

Specified Nodal Degree of Freedom Card(s) 


Column 


Contents 


“ITT” 

ND 

Node Number of Specified D.O.F, 


6-10 

IDR 

Direction of Specified D.O.F. 




(1 = u, 2 = V, and 3 = w) 


11-20 

VBDF(N) 

Specified Displacement Value 



^Repeat Card 12 NSDF Times 



* If NSBF 

.EQ. 0 - Skip Card 13 and Go on to Card 14 

Card 

13 (215, DIO. 5) 

Specified Nodal Force Card(s) 


Column 


Contents 


1-5 

ND 

Node Number of Specified Force 


6-10 

IDR 

Direction of Specified Force 




(1 = u, 2 = v, 3 = w) 


11-20 

VBSF(N) 

Specified Boundary Force Value 



^Repeat Card 13 NSBF Times 

Card 

14 (8D10.5) 

Normal Strain and Hygrothermal Load Card 


Column 


Contents 


1-10 

EPSX 

Applied Normal, (X-Di recti on ) , Strain 


11-20 

TEMP 

Temperature Change or Moisture Concentration 


14 


5.0 A Note on Symmetry Requirements 

To reduce the cost of a finite element stress analysis symmetry 
arguments are often used when building the finite element mesh. To 
determine the response of a thin composite laminate tensile specimen the 
analyst need only model one quarter of the specimen cross-section. To 
obtain valid results, however, the appropriate symmetry boundary condi- 
tions must be specified. These boundary conditions are shown in 
Fig. (1.4). Along the z-axis (y=0) the U and V nodal displacements must 
be specified to be zero. Similarly, along the y-axis (z=0) the W dis- 
placements must be set equal to zero. 




ORIGINAL PAGE !9 
OF POOR QUALITY 


The upper right quadrant of the specimen cross-section will 
be modeled using CLFE2D 


r 

',z) = 0_ 



The appropriate boundary conditions for the quarter 
symmetry model are shown above 


Fig. 1.4 A ppropriate Boundary Conditions for a Quarter 
Symmetry Model 
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6.0 Example Problem 

Consider an eight ply T300/5208 Graphite Epoxy laminate with a 
[+45 2 /-452 ] s stacking sequence. The laminate is 0,50 inches wide and 
each ply is 0.005 inches thick. It is desired to study the response of 
the laminate to an applied tensile normal strain e° = 0.001. Of par- 
ticular interest are the stress, strain and strain energy density dis- 
tributions along the +45°/-45° interface. 

Solution : A 36 element, 42 node quarter symmetry mesh was developed to 

study the given problem using CLFE2D. The input file and output for 
this problem are shown in Appendix A. Note that the stress strain 
output code varies for each element in the mesh. By doing this, the 
overall response of the laminate and the response at the +45°, -45° 
interface can be obtained with a minimum of output. 

A simple finite element mesh was used in the solution of the 
example problem for the purpose of exposition. However, CLFE2D is 
capable of solving large scale research problems. 


o 
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origwal page is 
OF POOR QUALfTY 


EXAMPLE INPUT DATA 


EXAMPLE PROBLEM, T300/5208, 8-LAYER, ( 45 . 0 , 45 . 0, -45 , 0, -45 . 0 ) S 
EPSX= 0.001, HYGROTHERMAL LOAD = 0.000 


36 

42 

4 

16 

0 

3 

0 

0 

0 0 



0.01 

2 , 50D-03 







1 

1 

1 

1 

1 

1 






32.0 

100.0 


0.5 


0.5 



920D+08 

1 , 560D+06 

1 . 560D+06 

. 4870D+06 

. 8200D+06 

. 8200D+06 

0 

‘.490 

0 

.238 

0 

.238 


0.00 

. 1380D-04 

. 1380D-04 


45.0 


45.0 

- 

45.0 

- 

45.0 



1 

4 

9 

10 

5 

1 

4 




2 

3 

8 

9 

4 

2 

4 

6 

1 7 


3 

2 

7 

8 

3 

3 

4 

9 

3 8 


4 

1 

6 

7 

2 

4 

4 




5 

9 

14 

15 

10 

1 





6 

8 

13 

14 

9 

2 

6 

1 

7 


' 7 

7 

12 

13 

8 

3 

9 

3 

8 


8 

6 

11 

12 

7 

4 





9 

15 

23 

24 

15 

1 

5 




10 

14 

23 

15 

14 

1 

5 




11 

14 

22 

23 

14 

1 

5 




12 

14 

21 

22 

14 

2 

5 




13 

13 

21 

14 

13 

2 

5 




14 

13 

20 

21 

13 

2 

5 




15 

13 

19 

20 

13 

3 

5 




16 

12 

19 

13 

12 

3 

5 




17 

12 

18 

19 

12 

3 

5 




18 

12 

17 

18 

12 

4 

5 




19 

11 

17 

12 

11 

4 

5 




20 

11 

16 

17 

11 

4 

5 




21 

23 

32 

33 

24 

1 





22 

22 

31 

32 

23 

1 





23 

21 

30 

31 

22 

2 





24 

20 

29 

30 

21 

2 

6 

1 

7 


25 

19 

28 

29 

20 

3 

9 

3 

8 


26 

18 

27 

28 

19 

3 





27 

17 

26 

27 

18 

4 





28 

16 

25 

26 

17 

4 





29 

32 

41 

42 

33 

1 

5 

2 



30 

31 

40 

41 

32 

1 

5 




31 

30 

39 

40 

31 

2 

5 




32 

29 

38 

39 

30 

2 

5 

6 

1 7 


33 

28 

37 

38 

29 

3 

5 

9 

3 8 


3 4 

27 

3 6 

37 

28 

3 

5 




35 

26 

35 

36 

27 

4 

5 




36 

25 

34 

35 

26 

4 

5 




1 


0.0 


0.0 






2 


0.0 


2.0 






3 


0.0 


4.0 






4 


0.0 


6.0 






5 


0.0 


8.0 






6 


8.0 


0.0 







■Vfc. * ^ 


20 


7 


8.0 


2.0 

8 


8.0 


4.0 

9 


8.0 


6.0 

10 


8.0 


8.0 

11 


16.0 


0.0 

12 


16.0 


2.0 

13 


16.0 


4.0 

14 


16.0 


6.0 

15 


16.0 


8.0 

16 


20.0 


0.0 

17 


20.0 


1.0 

18 


20.0 


2.0 

19 


20.0 


3.0 

20 


20.0 


4.0 

21 


20.0 


5.0 

22 


2 G .0 


6,0 

23 


2 C . 0 


7.0 

24 


20.0 


8.0 

25 


23.0 


0.0 

26 


23.0 


1.0 

27 


23.0 


2.0 

28 


23.0 


3.0 

29 


23.0 


4.0 

30 


23.0 


5.0 

31 


23.0 


6.0 

32 


23.0 


7.0 

33 


23.0 


8.0 

34 


25.0 


0.0 

35 


25.0 


1.0 

36 


25 . 0 


2.0 

37 


25.0 


3.0 

38 


25.0 


4.0 

39 


25.0 


5.0 

40 


25.0 


6.0 

41 - 


25.0 


7.0 

42 


25.0 


8.0 

1 

1 


0.0 


2 

1 


0.0 


3 

1 


0.0 


4 

1 


0.0 


5 

1 

* 

0.0 


1 

2 


0.0 


2 

2 


0.0 


3 

2 


0.0 


4 

2 


0.0 


5 

2 


0.0 


1 

3 


0.0 


6 

3 


0.0 


11 

3 


0.0 


16 

3 


0.0 


25 

3 


0.0 


34 

3 


0.0 




0 

.001 





0 . DO 
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Example Problem, T300/5208, 8-Layer, (45.0,45.0 
Epsx = 0.001, Hygrothermal Load = 0.000 


